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ABSTRACT 

We report low- and high- resolution spectra of comet C/2002C1 (Ikeya-Zhang) from McDonald Obser- 
vatory. The comet had a well-developed ion tail including CO + , CO^, CH+, and H 2 + . We used our 
high-resolution spectra to search for Nj. None was detected and we placed upper limits on N^/CO + of 
5.4 x 1CP 4 . N^ was detected in the low-resolution spectra but we show that this emission was probably 
telluric in origin (if cometary, we derive N^/CO + = 5.5 x 10~ 3 , still very low). We discuss the implica- 
tions for the conditions in the early solar nebula of the non-detection of N J . These depend on whether 
the H2O ice was deposited in the amorphous or crystalline form. If H2O was deposited in its crystalline 
form, the detection of CO + but not N^ has implications for H2O/H2 in the early solar nebula. 

Subject headings: Comets: individual(Ikeya-Zhang) — Solar System: formation 
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1. INTRODUCTION 

Knowledge of the nitrogen content of comets is impor- 
tant for an understanding of conditions in the early solar 
nebula. Nitrogen exists in cometary ices in many forms 
including N2, NH3, HCN, etc. Most of these species are 
chemically reactive in the coma gases making it difficult to 
use coma observations to unravel the nitrogen chemistry 
of the solar nebula. However, conditions in the early solar 
nebula were such that the dominant equilibrium species of 
carbon, oxygen and nitrogen should be N 2 , CO, and H2O 
(Lewis & Prinn 1980). Observations of N2H + in dense 
molecular clouds, coupled with chemical models, led Wom- 
ack et al. (1992) to infer that in potential star-forming 
regions nitrogen is preferentially in N2, rather than NH3. 

N2 is the least reactive of the nitrogen-bearing species 
and is thus the most appropriate to study to understand 
the nitrogen chemistry of comets. However, observa- 
tions of cometary N 2 are exceedingly difficult. Ground- 
based spectra are hampered by the telluric N 2 atmosphere; 
spacecraft flyby mass spectrometer observations are com- 
promised by the fact that N2 shares the mass 28 bin with 
CO, which is known to be quite common in cometary co- 
mae. 

A suitable method for measuring the N2 content of 
cometary ices is by studying its ion, N^~. Generally, ob- 
servations are obtained of the Nj IN B 2 S+ - A : S+ (0,0) 
band at 3914A. Successful observations of this band re- 
quire that the comet have a well-developed ion tail and 
that the data be obtained with sufficiently high spectral 
resolution to isolate the band from other cometary emis- 
sions and from any telluric N^ emission. 

Comet C/2002C1 (Ikeya-Zhang) was discovered in early 
February 2002 and reached perihelion on 2002 March 18 
at a heliocentric distance of 0.507AU. Our spectral ob- 
servations showed that it had a strong ion tail containing 
ions of CO+, COj, CH+, and H 2 0+. The comet was 
relatively bright, making it possible to observe at high 
spectral-resolution. Thus, we obtained spectra of Ikeya- 



Zhang at both high- and low-resolution in order to search 
for the signature of Nj in the tail. In this paper, we report 
on our non-detection of any Nj attributable to the comet 
and discuss the implications of our derived upper limits. 

2. OBSERVATIONS 

We obtained spectra of comet Ikeya-Zhang during April 
2002 using the 2.7-m Harlan J. Smith Telescope at Mc- 
Donald Observatory. Table 1 gives the parameters of 
the observations. Two instruments were used: the Large 
Cassegrain Spectrograph (LCS) is a long-slit CCD in- 
strument with a resolving power of R=A/AA=550; the 
2dcoude is a cross-dispersed echelle spectrograph with 
R=60,000. The LCS observations covered the spectral 
bandpass from 3000-5600A. The slit was 2 arcsec wide and 
approximately 2.5arcmin long, with each pixel imaging 
1.28 arcsec on the sky. The 2dcoude observations covered 
the spectral region from 3700-5700A continuously, with 
coverage to 1.02/im with increasing interorder gaps. The 
slit was 1.2 x 8.2 arcsec. With both instruments, observa- 
tions were obtained of the optocenter, as well as of the tail. 
The LCS data reduction procedure has been described in 
detail by Cochran et al. (1992); the 2dcoude data reduc- 
tion was described in Cochran et al. (2000). 

Figure 1 shows spectra obtained with the LCS, both on 
the optocenter and 31,400 km tailward. The ion spectrum 
of the tail is quite rich. Note that there is a feature to the 
red of the CN Av=0 band which is labeled as N^ in this 
figure. Is this really N^ emission? Some of it is certainly 
contamination from the CO + (5,1) band. However, the 
relative fluorescence efficiencies in Table 3 of Magnani & 
A'Hearn (1986) can be used to predict that CO + from the 
(5,1) band contributes no more than 1/4 of the detected 
feature. Another possible contaminant is an unidentified 
band of CO^ with bandheads at 3906 and 3916A. Evi- 
dence from Mrozowski (1947a, b) suggests that the con- 
tamination from this band is even less than that of the 
CO + (5,1) band. Therefore, most of the observed feature 
at 3914A can be attributed to N J . Is this N^ produced in 
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the cometary coma or is it telluric in origin? We measured 
the integrated band fluxes of this feature, along with the 
ionic CO + (3,0) and the neutral CN Av=-1 bands as a 
function of position within the coma. While the CN and 
CO + features are fairly isolated, making them straight- 
forward to measure, the possible feature sits on the 
wing of the very strong CN Av=0 band. Thus, extreme 
care must be taken when removing the underlying "contin- 
uum" of the feature. This was done by fitting the spectra 
at each position interactively. 

Figure 2 shows the derived integrated band fluxes as a 
function of position within the coma. The CN parent un- 
dergoes photodissociation as it flows outwards from the 
nucleus. In turn, the CN is photodissociated so that there 
is a monotonic decrease with cometocentric distance for 
the neutral CN. The CO + shows an increasing flux until 
about 25,000 km from the optocenter and then a slight de- 
cline. This is a result of ionization within the coma and 
entrainment in the solar wind. The possible N^" feature is 
mostly constant across the slit, unlike either the neutral 
CN or the ionic CO + . The S/N is low for this feature and 
removing the underlying CN adds uncertainty. If, as we 
believe, the feature is indeed constant across the slit, it 
would be the result of the terrestrial atmosphere. How- 
ever, we cannot rule out that N^ increases similarly to the 
CO+, implying that there is a cometary contributor. 

High-resolution spectra should be more sensitive than 
low-resolution spectra to weak line emission. Therefore, 
on 2002 April 22, we obtained echelle spectra to search 
for the N^ (0,0) band. Figure 3 shows our tail spectrum 
of this spectral region. The expected location of N^ (0,0) 
band emission lines are marked above the spectrum. In- 
spection of the figure reveals no obvious line emission at 
the wavelengths marked. There is, however, an obvious 
broad feature at around 3914. 5A. This feature begins just 
to the red of the Nj P-branch bandhead. A similar, but 
weaker, feature was seen in comparable spectra of comet 
122P/deVico (Cochran et al. 2000). In that paper we ar- 
gued that the feature was the neutral C3 (0,2,0)-(0,0,0) 
band. In the deVico spectra, there was a single feature, 
while close inspection of Figure 3 reveals that in the Ikeya- 
Zhang spectrum we have actually detected an additional 
feature at around 3916A. The first feature coincides with 
the C3 R-branch bandhead while the second coincides with 
theQ-branch. The presence of two of the C 3 (0,2,0)-(0,0,0) 
bandheads confirms our identification of these features as 
C3. Indeed, comparison with the optocenter 2dcoude spec- 
tra taken both before and after the tail spectrum confirm 
this identification as well. 

3. UPPER LIMITS 

A limitation of the high-resolution 2dcoude data is that 
it is not easy to calibrate absolutely the spectra. However, 
we can calibrate the relative flux of different orders so we 
can derive ratios of abundances of species for comparison 
with models. To compute the upper limits of Nj /CO + , we 
need to measure an upper limit for the Nj non-detection 
and the integrated intensity of the CO + (2,0) band. 

For the integrated CO + intensity, we fit a continuum to 
the region of the detected lines and summed the counts in 
the band above the continuum. Then, the contribution of 
the strong CH + P(3) line was removed. Finally, we mul- 



tiplied our integrated band intensity by two to reflect the 
fact that we were observing only the 2 n 1 / 2 (F2) branches. 
The integrated band intensity was 634 counts. 

The N^ upper limit is determined by asking how much 
signal would we be able to hide within the noise. We mea- 
sured the rms in the bandpass of the P-branch. Then, the 
upper limit is just 1/2 x rms x bandpass, in appropriate 
units. These are 2a upper limits. To account for the com- 
plete Nj band, we multiplied the P-branch upper limits 
by 2, assuming a comparable R-branch upper limit. This 
results in a 2a upper limit of 4.2 counts. 

To convert from the two band intensities to a ratio of 
the abundances, we need to account for the relative re- 
sponse of the instrument at the two bandpasses. Using 
solar spectra obtained with the same instrument, we de- 
termined that we needed to multiply the counts in the 
order containing N^ by a factor of 1.6. In addition, 
we need to account for the relative fluorescence efficien- 
cies of the two species. We used excitation factors of 
7.0 xl0~ 2 photons sec -1 moU 1 for the N^ (0,0) band (Lutz 
et al. 1993) and 3.55 x 10~ 3 photons sec -1 moU 1 for the 
CO + (2,0) band (the average value from Figure 2 of Mag- 
nani & A'Hearn (1986)). Applying these factors, we find 
that 



comparable to upper limits found for deVico and Hale- 
Bopp with the same instrument (Cochran et al. 2000). 

Using the flux of the Nj feature detected in the low- 
resolution spectra, we can derive a ratio of N^ /CO + = 
5.5 x 10~ 3 , a factor of 10 greater than our upper limit 
derived from the high-resolution spectra. This demon- 
strates that we would have easily detected the signal seen 
in the low-resolution spectra had it been present in the 
high-resolution spectra, regardless of whether it was tel- 
luric or cometary. If the Nj in the low-resolution spectra 
is cometary, it indicates a very variable parent. 

If the N^ in the low-resolution spectra was telluric, why 
didn't we see it in the high-resolution spectra at the appro- 
priate Dopplcr shift? The N^ B 2 T<+ state requires 18.6eV 
to be populated. Therefore, in the Earth's atmosphere, it 
can only be found in dayglow or auroral spectra. There are 
no photochemical reactions which can produce this emis- 
sion in nightglow (Torr et al. 1992). The geomagnetic K p 
indices for the nights of our observations were very simi- 
lar and showed relatively quiescent ionospheric conditions. 
Thus, the presence of the telluric N^ in the low-resolution 
data cannot be the result of auroral activity. 

The explanation for the presence of telluric N^ in only 
some of the spectra lies in the position of the comet in 
the sky at the times of the observations. While both the 
low- and high-resolution spectra were obtained at approx- 
imately the same hour of the night, around astronomical 
twilight, the orbit of the comet results in our observing 
the comet at very different zenith distances during the 
two observing runs. During the early April timeframe of 
the low-resolution data, the zenith distance of the comet 
was between 70 and 80 degrees; during the 2dcoudc ob- 
servations, the zenith distance was between 42 and 52 de- 
grees. Recall that the telluric Nj emission is a dayglow 
feature. It is also excited high in the ionosphere, with only 
5% of the N^ vibrationally excited at 100 km and 50% at 
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450km (Fox & Dalgarno 1985). Thus, we could only de- 
tect it when the slant-path of the Sun illuminated the part 
of the atmosphere we were observing through, or only at 
very high zenith distances. Thus, if the was telluric, it 
is consistent that we only obseved it in the low-resolution, 
high-zenith distance observations. 

4. IMPLICATIONS 

Ikeya-Zhang represents the third comet for which we 
have obtained very high-quality, high-resolution spectra 
in which we did not detect and could place tight con- 
straints on the relative abundances of N^ and CO + . Even 
if our low resolution detections were of cometary ori- 
gin, we have observed very little N^ relative to CO + and 
this ratio is variable for a single comet. Previous detec- 
tions and upper-limits from other observers were discussed 
in Cochran et al. (2000). In general, almost none of these 
previous observations were obtained at comparable reso- 
lutions and high sensitivity as our observations. 

There is an expectation that neutral N2 and CO should 
be very common in the early solar nebula (Lewis & Prinn 
1980) and thus, they are quite likely deposited in ices in 
proportion to their source. The quantity of these species 
deposited in the ices is dependent on several factors, in- 
cluding the deposition temperature of the ice. Current 
models of the dynamical evolution of the early Solar Sys- 
tem suggest that the Oort cloud comets we have studied 
were formed in the Uranus-Neptune region prior to pertur- 
bation into the Oort cloud (c.f. Weissman (1991); Duncan 
et al. (1987)). The temperature in this region was probably 
around 50±20K (Boss et al. 1989). Indeed, evidence con- 
firms this deposition temperature for Oort cloud comets 
(Meier et al. 1998a,b). 

The implications for the expected quantities of CO and 
N2 ices in cometary nuclei is then dependent on what form 
of H 2 ice is deposited in the nucleus. If the solar neb- 
ula ices never experienced warm temperatures, the H 2 
ice would be amorphous at the time of the deposition of 
the nucleus. Laboratory experiments of the deposition of 
CO and N 2 gas in amorphous H 2 ices have shown that 
CO is trapped 15-30 times more efficiently than N 2 in 
amorphous ice if the temperature of the amorphous ice 
is around the temperatures expected for the solar nebula 
(Bar-Nun et al. 1988; Notesco & Bar-Nun 1996). Using 
these results, Owen & Bar-Nun (1995) concluded that icy 
planctcsimals which formed in the Uranus-Neptune region 



would have N 2 /CO« 0.06 in the gases trapped in the ice if 
N 2 /CO« 1 in the nebula. This value is much higher than 
any of our upper limits (or pssible detection) for comets 
Ikeya-Zhang, deVico or Hale-Bopp. 

Mousis et al. (2000) suggested that the H 2 was not 
deposited in its amorphous state because of turbulence in 
the nebula at the time of deposition. Instead, the water ice 
infalling in the presolar cloud would first have vaporized 
(Chick & Cassen 1997) and subsequently been condensed 
in the nebula in the form of crystalline ice (Kouchi et al. 
1994). As the nebula continued to cool and H 2 ice crys- 
tallized, the other volatiles would be trapped in hydrate 
clathrates (Iro et al. 2002). 

The trapping of the other volatiles via clathration is de- 
pendent on the quantity of H 2 ice since clathration con- 
sumes water. Iro et al. (2002) pointed out that if H 2 0/H 2 
when water condensed was less than 2.8 times the solar 
O/H ratio, the clathration process would stop prior to 
trapping all of the volatiles. Then, the ratio of the trapped 
volatiles would be dependent on the curve of stability for 
that species (Lunine & Stevenson 1985). Inspection of Fig- 
ure 2 of Iro et al. (2002) shows that CO would condense 
at warmer temperatures than N 2 . Thus, Iro et al. (2002) 
suggested that our upper limits for N 2 f /CO + imply that 
H 2 0/H 2 < 2.8 x solar O/H when deVico and Hale-Bopp 
condensed. Our new upper limit for N+/CO+ obtained 
for Ikeya-Zhang is consistent with this picture. 

The three comets for which we have firm upper limits for 
N^" / CO + are consistent with formation of cometary nuclei 
from crystalline H 2 ice and subsequent incomplete trap- 
ping of volatiles in hydrate clathrates. All three of the 
comets are Oort cloud comets. Iro et al. (2002) pointed 
out that conditions for comets formed further from the 
proto-Sun, e.g., in the Kuiper belt, might allow for the 
deposition of amorphous H 2 ice. Thus, it would be de- 
sirable to obtain high-resolution spectra of the ion tail of 
any Kuiper belt comet which is bright enough in the future 
to test our understanding of the conditions for deposition 
of ices in the outer solar nebula. However, such observa- 
tions might not be conclusive because short-period comets 
have undergone a great deal of volatile-depleting process- 
ing during repeated solar passages. 

This work was supported by NASA Grant NAG5-9003. 
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Table 1 
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Fig. 1. — Low-resolution spectra of comet Ikeya-Zhang. Shown are spectra of the optocenter and 31,400 km tailward of the optocenter. In 
the upper panel, both spectra are shown to the same scale but offset from one another (the scale for the optocenter is on the left; that for the 
tail is on the right). Many features are marked. In the lower panel the ordinate of the tail spectrum has been expanded to show the weak 
ionic features. The text discusses the feature marked N^". 
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Fig. 2. — The distribution of the gas within the coma. Shown are examples of the distribution of a neutral (CN) and an ionized (CO + ) 
species as a function of cometocentric distance. In addition, the integrated band flux of the feature denoted N^" in Figure 1 is also shown. 
The distribution of this species is flat, unlike the CN or CO+ gas. The flat distribution is a signature of a telluric emission feature. 
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Fig. 3. — High-resolution spectra of comet Ikeya-Zhang. Shown is the spectrum of the tail in the region of the Nj (0,0) band. The expected 
positions of the P- and R-branch lines are marked. None are detected. Also marked are the expected positions of the C3 (0,2,0)-(0,0,0) 
band emissions. The R- and Q-branches of C3 were detected, even 15,000 km from the optocenter. 



